Rare CRHR2 and GRM8 variants identified as candidate factors associated with eating disorders in Japanese patients by whole exome sequencing.

Eating disorders (EDs) are a type of psychiatric disorder characterized by pathological eating and related behavior and considered to be highly heritable. The purpose of this study was to explore rare variants expected to display biological functions associated with the etiology of EDs. We performed whole exome sequencing (WES) of affected sib-pairs corresponding to disease subtype through their lifetime and their parents. From those results, rare single nucleotide variants (SNVs) concordant with sib-pairs were extracted and estimated to be most deleterious in the examined families. Two non-synonymous SNVs located on corticotropin-releasing hormone receptor 2 (CRHR2) and glutamate metabotropic receptor 8 (GRM8) were identified as candidate disease susceptibility factors. The SNV of CRHR2 was included within the cholesterol binding motif of the transmembrane helix region, while the SNV of GRM8 was found to contribute to hydrogen bonds for an α-helix structure. CRHR2 plays important roles in the serotoninergic system of dorsal raphe nuclei, which is involved with feeding and stress-coping behavior, whereas GRM8 modulates glutamatergic neurotransmission. Moreover, GRM8 modulates glutamatergic neurotransmission, and is also considered to have effects on dopaminergic and adrenergic neurotransmission. Thus, identification of rare and deleterious variants in this study is expected to increase understanding and treatment of affected individuals. Further investigation regarding the biological function of these variants may provide an opportunity to elucidate the pathogenesis of EDs.

Two-step evolution of HIV-1 budding system leading to pandemic in the human population.

The pandemic HIV-1, HIV-1 group M, emerged from a single spillover event of its ancestral lentivirus from a chimpanzee. During human-to-human spread worldwide, HIV-1 diversified into multiple subtypes. Here, our interdisciplinary investigation mainly sheds light on the evolutionary scenario of the viral budding system of HIV-1 subtype C (HIV-1C), a most successfully spread subtype. Of the two amino acid motifs for HIV-1 budding, the P(T/S)AP and YPxL motifs, HIV-1C loses the YPxL motif. Our data imply that HIV-1C might lose this motif to evade immune pressure. Additionally, the P(T/S)AP motif is duplicated dependently of the level of HIV-1 spread in the human population, and >20% of HIV-1C harbored the duplicated P(T/S)AP motif. We further show that the duplication of the P(T/S)AP motif is caused by the expansion of the CTG triplet repeat. Altogether, our results suggest that HIV-1 has experienced a two-step evolution of the viral budding process during human-to-human spread worldwide.

Genomic diversity of SARS-CoV-2 can be accelerated by mutations in the nsp14 gene.

Coronaviruses, including severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), encode a proofreading exonuclease, nonstructural protein 14 (nsp14), that helps ensure replication competence at a low evolutionary rate compared with other RNA viruses. In the current pandemic, SARS-CoV-2 has accumulated diverse genomic mutations including in nsp14. Here, to clarify whether amino acid substitutions in nsp14 affect the genomic diversity and evolution of SARS-CoV-2, we searched for amino acid substitutions in nature that may interfere with nsp14 function. We found that viruses carrying a proline-to-leucine change at position 203 (P203L) have a high evolutionary rate and that a recombinant SARS-CoV-2 virus with the P203L mutation acquired more diverse genomic mutations than wild-type virus during its replication in hamsters. Our findings suggest that substitutions, such as P203L, in nsp14 may accelerate the genomic diversity of SARS-CoV-2, contributing to virus evolution during the pandemic.

Comprehensive Investigation on the Interplay between Feline APOBEC3Z3 Proteins and Feline Immunodeficiency Virus Vif Proteins.

As the hosts of lentiviruses, almost 40 species of felids (family Felidae) are distributed around the world, and more than 20 feline species test positive for feline immunodeficiency virus (FIV), a lineage of lentiviruses. These observations suggest that FIVs globally infected a variety of feline species through multiple cross-species transmission events during a million-year history. Cellular restriction factors potentially inhibit lentiviral replication and limit cross-species lentiviral transmission, and cellular APOBEC3 deaminases are known as a potent restriction factor. In contrast, lentiviruses have evolutionary-acquired viral infectivity factor (Vif) to neutralize the APOBEC3-mediated antiviral effect. Because the APOBEC3-Vif interaction is strictly specific for viruses and their hosts, a comprehensive investigation focusing on Vif-APOBEC3 interplay can provide clues that will elucidate the roles of this virus-host interplay on cross-species transmission of lentiviruses. Here, we performed a comprehensive investigation with 144 patterns of a round robin test using 18 feline APOBEC3Z3 genes, an antiviral APOBEC3 gene in felid, and 8 FIV Vifs and derived a matrix showing the interplay between feline APOBEC3Z3 and FIV Vif. We particularly focused on the interplay between the APOBEC3Z3 of three felids (domestic cat, ocelot, and Asian golden cat) and an FIV Vif (strain Petaluma), and revealed that residues 65 and 66 of the APOBEC3Z3 protein of multiple felids are responsible for the counteraction triggered by FIV Petaluma Vif. Altogether, our findings can be a clue to elucidate not only the scenarios of the cross-species transmissions of FIVs in felids but also the evolutionary interaction between mammals and lentiviruses. IMPORTANCE Most of the emergences of new virus infections originate from the cross-species transmission of viruses. The fact that some virus infections are strictly specific for the host species indicates that certain "species barriers" in the hosts restrict cross-species jump of viruses, while viruses have evolutionary acquired their own "arms" to overcome/antagonize/neutralize these hurdles. Therefore, understanding of the molecular mechanism leading to successful cross-species viral transmission is crucial for considering the menus of the emergence of novel pathogenic viruses. In the field of retrovirology, APOBEC3-Vif interaction is a well-studied example of the battles between hosts and viruses. Here, we determined the sequences of 11 novel feline APOBEC3Z3 genes and demonstrated that all 18 different feline APOBEC3Z3 proteins tested exhibit anti-feline immunodeficiency virus (FIV) activity. Our comprehensive investigation focusing on the interplay between feline APOBEC3 and FIV Vif can be a clue to elucidate the scenarios of the cross-species transmissions of FIVs in felids.

TRIM28 Expression on Dendritic Cells Prevents Excessive T Cell Priming by Silencing Endogenous Retrovirus.

Acquired immune reaction is initiated by dendritic cells (DCs), which present Ags to a few naive Ag-specific T cells. Deregulation of gene expression in DCs may alter the outcome of the immune response toward immunodeficiency and/or autoimmune diseases. Expression of TRIM28, a nuclear protein that mediates gene silencing through heterochromatin, decreased in DCs from old mice, suggesting alteration of gene regulation. Mice specifically lacking TRIM28 in DCs show increased DC population in the spleen and enhanced T cell priming toward inflammatory effector T cells, leading to acceleration and exacerbation in experimental autoimmune encephalomyelitis. TRIM28-deficient DCs were found to ectopically transcribe endogenous retrovirus (ERV) elements. Combined genome-wide analysis revealed a strong colocalization among the decreased repressive histone mark H3K9me3-transcribed ERV elements and the derepressed host genes that were related to inflammation in TRIM28-deficient DCs. This suggests that TRIM28 occupancy of ERV elements critically represses expression of proximal inflammatory genes on the genome. We propose that gene silencing through repressive histone modification by TRIM28 plays a role in maintaining the integrity of precise gene regulation in DCs, which prevents aberrant T cell priming to inflammatory effector T cells.

Comprehensive genomic analysis reveals dynamic evolution of endogenous retroviruses that code for retroviral-like protein domains.

BACKGROUND: Endogenous retroviruses (ERVs) are remnants of ancient retroviral infections of mammalian germline cells. A large proportion of ERVs lose their open reading frames (ORFs), while others retain them and become exapted by the host species. However, it remains unclear what proportion of ERVs possess ORFs (ERV-ORFs), become transcribed, and serve as candidates for co-opted genes. RESULTS: We investigated characteristics of 176,401 ERV-ORFs containing retroviral-like protein domains (gag, pro, pol, and env) in 19 mammalian genomes. The fractions of ERVs possessing ORFs were overall small (~ 0.15%) although they varied depending on domain types as well as species. The observed divergence of ERV-ORF from their consensus sequences showed bimodal distributions, suggesting that a large proportion of ERV-ORFs either recently, or anciently, inserted themselves into mammalian genomes. Alternatively, very few ERVs lacking ORFs were found to exhibit similar divergence patterns. To identify candidates for ERV-derived genes, we estimated the ratio of non-synonymous to synonymous substitution rates (dN/dS) for ERV-ORFs in human and non-human mammalian pairs, and found that approximately 42% of the ERV-ORFs showed dN/dS < 1. Further, using functional genomics data including transcriptome sequencing, we determined that approximately 9.7% of these selected ERV-ORFs exhibited transcriptional potential. CONCLUSIONS: These results suggest that purifying selection operates on a certain portion of ERV-ORFs, some of which may correspond to uncharacterized functional genes hidden within mammalian genomes. Together, our analyses suggest that more ERV-ORFs may be co-opted in a host-species specific manner than we currently know, which are likely to have contributed to mammalian evolution and diversification.

A role for gorilla APOBEC3G in shaping lentivirus evolution including transmission to humans.

The APOBEC3 deaminases are potent inhibitors of virus replication and barriers to cross-species transmission. For simian immunodeficiency virus (SIV) to transmit to a new primate host, as happened multiple times to seed the ongoing HIV-1 epidemic, the viral infectivity factor (Vif) must be capable of neutralizing the APOBEC3 enzymes of the new host. Although much is known about current interactions of HIV-1 Vif and human APOBEC3s, the evolutionary changes in SIV Vif required for transmission from chimpanzees to gorillas and ultimately to humans are poorly understood. Here, we demonstrate that gorilla APOBEC3G is a factor with the potential to hamper SIV transmission from chimpanzees to gorillas. Gain-of-function experiments using SIVcpzPtt Vif revealed that this barrier could be overcome by a single Vif acidic amino acid substitution (M16E). Moreover, degradation of gorilla APOBEC3F is induced by Vif through a mechanism that is distinct from that of human APOBEC3F. Thus, our findings identify virus adaptations in gorillas that preceded and may have facilitated transmission to humans.

Sequence Compression Benchmark (SCB) database-A comprehensive evaluation of reference-free compressors for FASTA-formatted sequences.

BACKGROUND: Nearly all molecular sequence databases currently use gzip for data compression. Ongoing rapid accumulation of stored data calls for a more efficient compression tool. Although numerous compressors exist, both specialized and general-purpose, choosing one of them was difficult because no comprehensive analysis of their comparative advantages for sequence compression was available. FINDINGS: We systematically benchmarked 430 settings of 48 compressors (including 29 specialized sequence compressors and 19 general-purpose compressors) on representative FASTA-formatted datasets of DNA, RNA, and protein sequences. Each compressor was evaluated on 17 performance measures, including compression strength, as well as time and memory required for compression and decompression. We used 27 test datasets including individual genomes of various sizes, DNA and RNA datasets, and standard protein datasets. We summarized the results as the Sequence Compression Benchmark database (SCB database, http://kirr.dyndns.org/sequence-compression-benchmark/), which allows custom visualizations to be built for selected subsets of benchmark results. CONCLUSION: We found that modern compressors offer a large improvement in compactness and speed compared to gzip. Our benchmark allows compressors and their settings to be compared using a variety of performance measures, offering the opportunity to select the optimal compressor on the basis of the data type and usage scenario specific to a particular application.

Alopecia areata susceptibility variant in MHC region impacts expressions of genes contributing to hair keratinization and is involved in hair loss.

BACKGROUND: Alopecia areata (AA) is considered a highly heritable, T-cell-mediated autoimmune disease of the hair follicle. However, no convincing susceptibility gene has yet been pinpointed in the major histocompatibility complex (MHC), a genome region known to be associated with AA as compared to other regions. METHODS: We engineered mice carrying AA risk allele identified by haplotype sequencing for the MHC region using allele-specific genome editing with the CRISPR/Cas9 system. Finally, we performed functional evaluations in the mice and AA patients with and without the risk allele. FINDINGS: We identified a variant (rs142986308, p.Arg587Trp) in the coiled-coil alpha-helical rod protein 1 (CCHCR1) gene as the only non-synonymous variant in the AA risk haplotype. Furthermore, mice engineered to carry the risk allele displayed a hair loss phenotype. Transcriptomics further identified CCHCR1 as a novel component interacting with hair cortex keratin in hair shafts. Both, these alopecic mice and AA patients with the risk allele displayed morphologically impaired hair and comparable differential expression of hair-related genes, including hair keratin and keratin-associated proteins (KRTAPs). INTERPRETATION: Our results implicate CCHCR1 with the risk allele in a previously unidentified subtype of AA based on aberrant keratinization in addition to autoimmune events. FUNDING: This work was supported by JSPS KAKENHI (JP16K10177) and the NIHR UCLH Biomedical Research center (BRC84/CN/SB/5984).

Efficient differentiation and polarization of primary cultured neurons on poly(lactic acid) scaffolds with microgrooved structures.

Synthetic biodegradable polymers including poly(lactic acid) (PLA) are attractive cell culture substrates because their surfaces can be micropatterned to support cell adhesion. The cell adhesion properties of a scaffold mainly depend on its surface chemical and structural features; however, it remains unclear how these characteristics affect the growth and differentiation of cultured cells or their gene expression. In this study, we fabricated two differently structured PLA nanosheets: flat and microgrooved. We assessed the growth and differentiation of mouse primary cultured cortical neurons on these two types of nanosheets after pre-coating with poly-D-lysine and vitronectin. Interestingly, prominent neurite bundles were formed along the grooves on the microgrooved nanosheets, whereas thin and randomly extended neurites were only observed on the flat nanosheets. Comparative RNA sequencing analyses revealed that the expression of genes related to postsynaptic density, dendritic shafts, and asymmetric synapses was significantly and consistently up-regulated in cells cultured on the microgrooved nanosheets when compared with those cultured on the flat nanosheets. These results indicate that microgrooved PLA nanosheets can provide a powerful means of establishing a culture system for the efficient and reproducible differentiation of neurons, which will facilitate future investigations of the molecular mechanisms underlying the pathogenesis of neurological disorders.

Nucleotide Archival Format (NAF) enables efficient lossless reference-free compression of DNA sequences.

SUMMARY: DNA sequence databases use compression such as gzip to reduce the required storage space and network transmission time. We describe Nucleotide Archival Format (NAF)-a new file format for lossless reference-free compression of FASTA and FASTQ-formatted nucleotide sequences. Nucleotide Archival Format compression ratio is comparable to the best DNA compressors, while providing dramatically faster decompression. We compared our format with DNA compressors: DELIMINATE and MFCompress, and with general purpose compressors: gzip, bzip2, xz, brotli and zstd. AVAILABILITY AND IMPLEMENTATION: NAF compressor and decompressor, as well as format specification are available at https://github.com/KirillKryukov/naf. Format specification is in public domain. Compressor and decompressor are open source under the zlib/libpng license, free for nearly any use. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Systematic survey of non-retroviral virus-like elements in eukaryotic genomes.

Endogenous viral elements (EVEs) are viral sequences that are endogenized in the host cell. Recently, several eukaryotic genomes have been shown to contain EVEs. To improve the understanding of EVEs in eukaryotes, we have developed a system for detecting EVE-like sequences in eukaryotes and conducted a large-scale nucleotide sequence similarity search using all available eukaryotic and viral genome assembly sequences (excluding those from retroviruses) stored in the National Center for Biotechnology Information genome database (as of August 14, 2017). We found that 3856 of 7007 viral genomes were similar to 4098 of 4102 eukaryotic genomes. For those EVE-like sequences, we constructed a database, Predicted Endogenous Viral Elements (pEVE, http://peve.med.u-tokai.ac.jp) which provides comprehensive search results summarized from an evolutionary viewpoint. A comparison of EVE-like sequences among closely related species may be useful to avoid false-positive hits. We believe that our search system and database will facilitate studies on EVEs.

Altered oligomeric states in pathogenic ALS2 variants associated with juvenile motor neuron diseases cause loss of ALS2-mediated endosomal function.

Familial amyotrophic lateral sclerosis type 2 (ALS2) is a juvenile autosomal recessive motor neuron disease caused by the mutations in the ALS2 gene. The ALS2 gene product, ALS2/alsin, forms a homophilic oligomer and acts as a guanine nucleotide-exchange factor (GEF) for the small GTPase Rab5. This oligomerization is crucial for both Rab5 activation and ALS2-mediated endosome fusion and maturation in cells. Recently, we have shown that pathogenic missense ALS2 mutants retaining the Rab5 GEF activity fail to properly localize to endosomes via Rac1-stimulated macropinocytosis. However, the molecular mechanisms underlying dysregulated distribution of ALS2 variants remain poorly understood. Therefore, we sought to clarify the relationship between intracellular localization and oligomeric states of pathogenic ALS2 variants. Upon Rac family small GTPase 1 (Rac1) activation, all mutants tested moved from the cytosol to membrane ruffles but not to macropinosomes and/or endosomes. Furthermore, most WT ALS2 complexes were tetramers. Importantly, the sizes of an ALS2 complex carrying missense mutations in the N terminus of the regulator of chromosome condensation 1-like domain (RLD) or in-frame deletion in the pleckstrin homology domain were shifted toward higher molecular weight, whereas the C-terminal vacuolar protein sorting 9 (VPS9) domain missense mutant existed as a smaller dimeric or trimeric smaller form. Furthermore, in silico mutagenesis analyses using the RLD protein structure in conjunction with a cycloheximide chase assay in vitro disclosed that these missense mutations led to a decrease in protein stability. Collectively, disorganized higher structures of ALS2 variants might explain their impaired endosomal localization and the stability, leading to loss of the ALS2 function.

New World feline APOBEC3 potently controls inter-genus lentiviral transmission.

BACKGROUND: The apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3; A3) gene family appears only in mammalian genomes. Some A3 proteins can be incorporated into progeny virions and inhibit lentiviral replication. In turn, the lentiviral viral infectivity factor (Vif) counteracts the A3-mediated antiviral effect by degrading A3 proteins. Recent investigations have suggested that lentiviral vif genes evolved to combat mammalian APOBEC3 proteins, and have further proposed that the Vif-A3 interaction may help determine the co-evolutionary history of cross-species lentiviral transmission in mammals. RESULTS: Here we address the co-evolutionary relationship between two New World felids, the puma (Puma concolor) and the bobcat (Lynx rufus), and their lentiviruses, which are designated puma lentiviruses (PLVs). We demonstrate that PLV-A Vif counteracts the antiviral action of APOBEC3Z3 (A3Z3) of both puma and bobcat, whereas PLV-B Vif counteracts only puma A3Z3. The species specificity of PLV-B Vif is irrespective of the phylogenic relationships of feline species in the genera Puma, Lynx and Acinonyx. We reveal that the amino acid at position 178 in the puma and bobcat A3Z3 is exposed on the protein surface and determines the sensitivity to PLV-B Vif-mediated degradation. Moreover, although both the puma and bobcat A3Z3 genes are polymorphic, their sensitivity/resistance to PLV Vif-mediated degradation is conserved. CONCLUSIONS: To the best of our knowledge, this is the first study suggesting that the host A3 protein potently controls inter-genus lentiviral transmission. Our findings provide the first evidence suggesting that the co-evolutionary arms race between lentiviruses and mammals has occurred in the New World.

Different effects of two mutations on the infectivity of Ebola virus glycoprotein in nine mammalian species.

Ebola virus (EBOV), which belongs to the genus Ebolavirus, causes a severe and often fatal infection in primates, including humans, whereas Reston virus (RESTV) only causes lethal disease in non-human primates. Two amino acids (aa) at positions 82 and 544 of the EBOV glycoprotein (GP) are involved in determining viral infectivity. However, it remains unclear how these two aa residues affect the infectivity of Ebolavirus species in various hosts. Here we performed viral pseudotyping experiments with EBOV and RESTV GP derivatives in 10 cell lines from 9 mammalian species. We demonstrated that isoleucine at position 544/545 increases viral infectivity in all host species, whereas valine at position 82/83 modulates viral infectivity, depending on the viral and host species. Structural modelling suggested that the former residue affects viral fusion, whereas the latter residue influences the interaction with the viral entry receptor, Niemann-Pick C1.

RXRB Is an MHC-Encoded Susceptibility Gene Associated with Anti-Topoisomerase I Antibody-Positive Systemic Sclerosis.

Systemic sclerosis is a systemic autoimmune and connective tissue disorder associated with the human leukocyte antigen locus. However, the functional relationship between human leukocyte antigen gene(s) and disease development remains unknown. To elucidate major histocompatibility complex-linked systemic sclerosis genetics, we performed genotyping of major histocompatibility complex-borne microsatellites and HLA-DPB1 alleles using DNA obtained from 318 anti-topoisomerase I antibody-positive patients and 561 healthy controls, all of Japanese descent. Those results revealed two major histocompatibility complex haplotypes associated with systemic sclerosis. Exome sequencing and targeted analysis of these risk haplotypes identified rs17847931 in RXRB as a susceptibility variant (P = 1.3 × 10-15; odds ratio [OR] = 9.4) with amino acid substitution p.V95A on the risk haplotype harboring HLA-DPB1∗13:01. No identical variant in the other haplotype including DPB1*09:01 was observed, though that haplotype also showed a significant association (P = 8.5 × 10-22; OR = 4.3) with systemic sclerosis. Furthermore, the number of risk factors was shown to be a predominant factor, as individuals with two factors had elevated risk (P = 6.7 × 10-13; OR = 30.2). We concluded that RXRB may be involved in antifibrotic activity in skin and chromatin remodeling.

Functional mutations in spike glycoprotein of Zaire ebolavirus associated with an increase in infection efficiency.

Ebola virus (EBOV) is extremely virulent, and its glycoprotein is necessary for viral entry. EBOV may adapt to its new host humans during outbreaks by acquiring mutations especially in glycoprotein, which allows EBOV to spread more efficiently. To identify these evolutionary selected mutations and examine their effects on viral infectivity, we used experimental-phylogenetic-structural interdisciplinary approaches. In evolutionary analysis of all available Zaire ebolavirus glycoprotein sequences, we detected two codon sites under positive selection, which are located near/within the region critical for the host-viral membrane fusion, namely alanine-to-valine and threonine-to-isoleucine mutations at 82 (A82V) and 544 (T544I), respectively. The fine-scale transmission dynamics of EBOV Makona variants that caused the 2014-2015 outbreak showed that A82V mutant was fixed in the population, whereas T544I was not. Furthermore, pseudotype assays for the Makona glycoprotein showed that the A82V mutation caused a small increase in viral infectivity compared with the T544I mutation. These findings suggest that mutation fixation in EBOV glycoprotein may be associated with their increased infectivity levels; the mutant with a moderate increase in infectivity will fix. Our findings showed that a driving force for Ebola virus evolution via glycoprotein may be a balance between costs and benefits of its virulence.